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gQ F 
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7 rpc 1 
P 

instantaneous buoyancy term 
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4c Tup 
P 

percentage by weight of pollutant material in the fuel from 



Table 1 
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molecular weight (g/mole) 
ambient pressure (mb) 

integral of the Gaussian probability function between minus 
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TK 

( z - z 
IK ml 

" t - 
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time after ignition required for the cloud to reach the stabiliza- 
tion height (sec) 


*R 


time after ignition (sec) 


si 
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adiabatic atmosphere (sec) 
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R Z mI 
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ml 

maximum rise of the ground cloud (sec) [ obtained from 
equation (l)] 

mean wind speed (m/sec) 
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height of stabilized cloud (m) 
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TECHNICAL MEMORANDUM X-3162 


DOWNWIND HAZARD CALCULATIONS FOR SPACE SHUTTLE 
LAUNCHES AT KENNEDY SPACE CENTER 
AND VANDENBERG AIR FORCE BASE 


SUMMARY 


The NA3A/MSFC Multilayer Diffusion Model is used to predict the 
dispersion from the rocket motors of the Space Shuttle at Kennedy Space Center 
and Vandenberg Air Force Base. The analysis of the dispersion of the rocket 
exhaust effluents Is for the nine meteorological regimes at Kennedy Space 
Center (KSC) and six at Vandenberg Air Force Base (VAFB). It is concluded 
that outside the Kennedy Space Center industrial area of 10 km and outside a 
similar distance at Vandenberg Air Force Base, if a 10-mln short term public 
limit of 4 ppm of HC1 is applied (NAS/NRC Committee on Toxicology), in none 
of the typical selected meteorological cases studied wat Sue 10-min limit 
exceeded. From related climatological studies meteorological conditions will 
exist at both launch sites that will constrain the shuttle launches. As an 
example, one case at KSC had a 5 ppm average concentration for a 10-min 
period at 2 km and 1. 9 at 10 km distance from the launch pad. Also, 
scavenging of the exhaust cloud by precipitation was not considered and may 
result in increased concentration* 


INTRODUCTION 
Purpose and Scope 

The purpose of this technical memorandum is to provide quantitative 
estimates of pollutant concentrations associated with the emission of three 
major combustion products, HC1, CO, and Al^, to the lower atmosphere dur- 
ing normal launches at Kennedy Space Center and Vandenberg Air Force Base. 
These estimates are provided by the Marshall Space Flight Center to assist the 
National Aeronautics and Space Administration in assessing the environmental 
Impact of Space Shuttle launch operations. Attention has been focused on cal- 
culations of pollutant concentrations near ground level at various downwind 
distances. Only calc'ilatlons for normal launches are presented in this report. 


* Anon. Potential Tropospheric Exhaust Cloud Constraint and Shuttle Launch 
Delay Risk Assessment: Internal Report by Aerospace Environment Division, 
Space Sciences Laboratory, NASA-Marshall Space Flight Center, August 1974. 


The pollutant concentration calculations have been made by using the 
computerized NASA/MSFC Multilayer Diffusion Models f 1, 2] In conjunction 
with appropriate emissions data for the Space Shuttle solid propellant engines 
and with meteorological data for a number of p elected situations characteristic 
of each launch site. In this report, results ai^ presented for nine sets of 
meteorological conditions at Kennedy Space Center and six sets at Vandenberg 
Air Force Base. 


Fuel Properties and Vehicle Rise Data 

Characteristic fuel properties used in the concentration calculations 
are given in Table 1 and are principally based on data from the ’'Mission 1 
Ascent Trajectories (MCR 200R1 Configuration)" dated June 18, 1973. 


Q 


F 


W • H + W • H 
s s L 


( 1 ) 


where 

(Qj = total heat output = (W_) 9.4512G x 10* • (H c ) 6.91 x io 2 

r oo 

+ (W ) 1.53116 X 10 6 • (H ) 5.00 X io 2 , 

L JLi 

Q_ = 7. 2964 X io 9 cal/sec 

r 


However, it is assumed that 1. 26 x 10 s kg/sec of water would be used 
to cool the launch pad and that all the water is vaporized by the engine. This 
results in a heat loss (Q ) of 7. 74 x io 8 cal/sec. The effective heat used in 

Li 

the plume rise calculation is 


Q E = Q F " Q L * 

(Q_) effective heat = (Q.,) 7.29640 x io 9 - (Q t ) 7.74 x IO 8 
ht r L 

Q e = 6. 52240 X 10* cal/sec 
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The liquid engine does not release any of the three pollutants, HC1, CO, 
and A^Oj. 

The altitude-time curve for the Space Shuttle is also required to cal- 
culate the rise of the ground cloud of exhaust products. A logarithmic least 
squares regression curve was fitted to the time height data for the first 40 sec 
of flight yielding the expression 


t D = 0.672341Z 0 * 48 284 , (3) 

K 

where t_ is time in seconds and Z is in meters. For convenience. Figure 1 
K 

shows a plot of vehicle altitude versus time calculated from equation ( 3) . 


Organization of the Report 

The following section contains a description of the cloud rise calcula- 
tions for the models. Cloud dimensions as vertical distributions of exhausts 
are discussed in the third section. A description of the MSFC Multilayer 
Diffusion Model used in the calculations of concentrations downwind from the 
launch is presented in the fourth section. The Inst section describes the 
results of the calculations at KSC and VAFB. 

Appendix A contains a derivation of the plume rise equations and 
Appendix B contains a description of the procedures used in specifying the 
source and meteorological inputs. Appendix C includes the input parameters 
for the NASA/MSFC Multilayer Diffusion Model, and the toxicity criteria are 
given in Appendix D. 


CLOUD RISE CALCULATIONS 
Normal Launches 

The burning of rocket engines during normal launches results in the 
formation of a cloud of hot exhaust products which subsequently rises and 
entrains ambient air until an equilibrium with ambient conditions is achieved. 
Experience in predicting the buoyant rise [3-14] from normal launches of 
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m 



A.._ - 


solid-fueled vehicles indicates the rise is best predicted using a cloud rise 
model for instantaneous sources. For solid-fueled vehicles residence times 
near the pad are relatively short. 

Each of the models for cloud height is subdivided into two categories 

to account for the atmospheric temperature lapse rate. The model assumes 

that the atmosphere is either quasi-adiabatic or stable. Here the quasi- 

adiabatic is where the adiabatic atmosphere is the. limit, which means that the 

potential temperature difference (A 0) is zero or less, where the potential 

temperature difference is given by A 0 = 0 , _ - 0 . .If 

max cloud height surface 

this potential temperature difference is positive, then the atmosphere is stable. 
Since in most cases of interest there will be an inversion layer present, the 
stable cloud rise formula is the normally utilized relation [ 15, 16] . 


Cloud Rise Models 

The maximum rise z^ for an instantaneous source as given by the 
expression' in a stable atmosphere is given by 



whereas, the maximum cloud rise z _ downwind from an instantaneous source 

ml 

in an adiabatic atmosphere is given by 


Z 


ml 



( 5 ) 


In deriving equation (4), it is assumed that the initial upward momentum 
imparted to the exhaust gases by reflection from the ground surface and launch 
pad hardware is insignificant in comparison with the effect of thermal buoyancy. 
Based on limited experience in predicting cloud rise from launches at 


a mm 


n 




4 

\ 



• P 
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Vandenberg Air Force Base, this assumption appears to be justified. The time 
required for the cloud to reach the stabilization height is given by the expression 


= 


1/2 


( 6 ) 


In calculating z _ from equation ( 5) , the instantaneous heat released 
ml 

Qj is obtained from the relationship 


Q I ' VR^ml 1 ’ 


( 7 ) 


An inspection of the equations given above reveals an interdependence 

between the calculated maximum cloud rise z T , the height over which the 

ml 

potential temperature gradient 9$/9z is measured, and the value of t„{z } 

K Dll 

used in obtaining . Thus, the final value of maximum cloud rise must be 

found through iteration of equations (4) or (5) . The height over which 94> /9 z 

is measured and the time t^{ z } are thus made consistent with the value of 

R m 

z _ calculated from the model, 
ml 

It should be noted that the height predictions for the cl cud rise assume 
uniform potential temperature gradients over the atmospheric region in which 
the equations are applied. If a discontinuous atmospheric environment exists, 
relative to that over the launch site, as the exhaust cloud moves downwind and 
develops, then significant variation in the potential temperature gradient (and 
cloud rise) may occur unless this is taken into consideration. This can occur 
in regions of geographic discontinuities such as the land-water interface at the 
Kennedy Space Center area. If adequate measurements of the atmospheric 
structure exist, then the appropriate inputs can be made into the cloud rise 
prediction equations. 
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CLOUD DIMENSIONS AND VERTICAL DISTRIBUTION OF 
EXHAUST PRODUCTS 


Source inputs required for the diffusion model calculations include the 
stabilization height of the exhaust cloud and cloud dimensions, ac well as the 
vertical distribution of exhaust products in the stabilized cloud. The calcula- 
tion of the stabilization height z m was described previously. The calculation 

of the dimensions of the stabilized cloud and the vertical distribution of exhaust 
products is described below. 


Dimensions of the Exhaust Cloud at Stabilization 

The general formula used to calculate the radius of the stabilized cloud 
at height z is given by the expression 


{ z> = 


yz 


; z s z 


m 


y (2z - z) > 200 m ; z ^ z 

' ' m m 


( 8 ) 


where 


y = (y r ) = 0. 64 


Note that for z > z , the minimum radius of the stabilized cloud is set equal 
to 200 m. m 

The cloud is assumed to be symmetrical about a vertical axis through 
the cloud centroid. The alongwind and crosswind source dimensions of the 
cloud in each of the layers are calculated under the following assumptions: 

• The distribution of exhaust products within the cloud is Gaussian 
in the plane of the horizon. 

• The concentration of exhaust products at a lateral distance of 1 
radius from the cloud vertical axis is 10 percent of the concentra- 
tion at the cloud axis. 
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The alongwind and crosswind source dimensions required for input to 
,..ie MSFC Diffusion Models are defined for each layer by 


"J K} 


a 



yz/2. 15 ; z ^ z 

ttl 

y (2z - z)/2. 15 & 93 m ; z2 z 

m m 


( 9 ) 


where 

z r = midpoint of the Kth layer 



The quantities z and are, respectively, the height of the top and base 

of the Kth \~j er. 

The corresponding vertical source dimension for each layer was cal- 
culated from the expression 



( 10 ) 


Equation (10) applies to a rectangular material distribution which has been 
assumed to apply along the vertical in the Kth layer. 


Calculation of the Vertical Source Strength Distribution 
in the Stabilized Exhaust Cloud 

The fraction of material by weight in each of the K layers, f{k) , for 
the launches was calculated from the expression 
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I QP{z TO > 


F{K> = 


Q( p {zTK^ " P ^ Z BK^ 


K = 1 
K > 1 


(ID 


P { z ) is the integral of the Gaussian probability function between minus 
BK 

infinity and the base of the Kth layer, z___ , and is equal to P( z DI . - z T /cr) . 

. , , f y _ _ _ BK BK ml 

or is equal to y[z = z^j}/2,15. 

The MSFC Diffusion Model described in the next section requires that 
source strength in each of the K layers be specified per unit height. Since 
the desired concentration units for HC1 and CO are parts per million, the com- 
plete expression for the source strength model input for the Kth layer is 



For Aip„ the desired concentration units are milligrams per cubic meter and 
the complete expression for source strength in the Kth layer is 


Q = V F W - t (Utm) 

K Z TK " Z BK \ g / 


(13) 


Equations (ll), (12), and (13) were used to obtain the model input values of 
O for the various meteorological regimes. 


NASA/MSFC MULTILAYER DIFFUSION MODEL 


The spatial description, in terms of concentration and dosage, of the 
dispersive transport of effluents from a discrete source is afforded by the 
NASA/MSFC Multilayer Diffusion Model. Specifically, this application of the 
model is for the prediction of toxicity distribution associated with the rocket 
exhaust effluents emitted during the launch of a space vehicle in order to assess 
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the resulting environmental impact. The dispersive description accorded by 
the Multilayer Diffusion Model is initiated at the point where the cloud of 
effluents reaches thermodynamic equilibrium with the environment, and there- 
fore depends strongly on the kinematic and thermodynamic profiles of the 
atmospheric conditions along with a knowledge of the exhaust effluents present 
in the cloud. 

The initial considerations in this section are given to the techniques of 
establishing the spatial location of the ground cloud equilibrium. Secondly, 
a general discussion of the Lagrangian dispersion of a point source is given. 
The final discussion in this section explains how the Multilayer Diffusion Model 
incorporates the general diffusion description to account for the two stages of 
exhaust emission and accounts for the environmental effects. The significant 
mathematical expressions supporting these discussions have been included in 
Appendices A, B, and C. 


Altitude of Cloud Equilibrium 

The effluent cloud rise relations are employed to determine the altitude 
at which the ground cloud reaches equilibrium with the environment. The 
importance of this location is that it serves as the origin of the dispersive 
description. This equilibrium point is chosen as the origin in order to eliminate 
complex thermodynamic considerations and to limit the problem solely to 
kinematics. 

The burning of rocket engines results in the formation of a cloud of hot 
exhaust products which subsequently rises and entrains ambient air until an 
equilibrium with ambient conditions is reached. For normal launches, this 
cloud is formed principally by the forced ascent of hot turbulent exhaust pro- 
ducts that have been deflected laterally and vertically by the launch pad hard- 
ware and the ground surface. The height at which this ground cloud stabilizes 
is determined by the vehicle type and atmospheric stability. The vehicle type 
determines whether a continuous or instantaneous source model is required. 

In the instantaneous source model, spherical entrainment is assumed; that is, 
the entrained ambient air enters the exhaust cloud uniformly from all directions. 
In the continuous source model, cylindrical entrainment is assumed; that is, the 
entrained ambient air enters the cloud uniformly only on the sides of the cylinder 
and not the ends. 


The generalized dosage model for a nearly instantaneous source is 
defined by the product of four modular terms: 

Dosage = { Peak Dosage Terms} x { Lateral Term} 
x {Vertical Term} x { Depletion Term} 


Thus, the mathematical description for the concentration and dosage 
models permit flexibility in application to various sources and for changing 
atmospheric parameters while always maintaining a rigorous mass balance. 

Two obvious differences exist. First, the peak concentration term 
refers to the concentration at the point x, y = 0, z = H (where x is the wind 
direction and H is any height) and is defined by the expression 


Q 


( 27 r) 3 / 2 <T x°y 


(14) 


where Q is the source strength and a is the standard deviation ot the con- 
centration distribution in the ith direction. The peak dosage term is given by 


Q 


27 ru cr a 

y z 


(15) 


where u is the moan wind speed. The second difference between these models 
is that the concentration contains a modular alongwind term to account for 
downstream temporal effects not considered in the dosage model. The along- 
wind term affords an exponential decay in concentration as a function of: cloud 
transit time, concentration distribution, and the mean wind speed. 

The lateral term, which is common to both models, is another exponen- 
tial decay term and is a function of the Gaussian spreading rate and the distance 
laterally from the mean wind azimuth. The vertical term, again common to 
both models, is a rather complex decay function since it contains a multiple 
reflection term for the point source which stops the vertical cloud development 
at the top of the mixing layer and eventually changes the form of the vertical 
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concentration distribution from Gaussian to rectangular. The last modular in 
both models is the depletion term. This term accounts for the loss of material 
by simple decay processes, precipitation scavenging, or gravitational settling. 


The Description of the Models in the NASA/MSFC 
Multilayer Diffusion Model 

The normal launch environment will usually involve an atmospheric 
structure comprised of several horizontal meteorological layers with distinctive 
wind velocity, temperature, and humidity regimes between the surface and a 5- 
km altitude. Large horizontal spatial variation in these meteorological param- 
eters may also occur in the surfac" layer as a consequence of changes in ter- 
rain or land-water interfaces, which is accounted for by the diffusion model. 

The general diffusion model for the concentration [ equation ( 14) ] and the 
dosage f equation ( 1 5) ] assumes an expanding volume about a moving point of 
reference in a homogeneous environment T 17 - 20] . 

To overcome the obvious shortcomings of the general diffusion model 
but to stay within the established bounds of classical fluid mechanics f 21] , a 
multiple layer concept is introduced to cope with the vertical and horizontal 
atmospheric gradients. Here, the general diffusion model Is applied to individ- 
ual horizontal layers in which the meteorological structure is reasonably 
homogeneous and independent of the neighboring layers. These layers have 
boundaries which are placed at points of major discontinuities in the vertical 
profiles of wind velocity, temperature, and humidity. Since the Multilayer 
Diffusion Model has imposed the general restriction of layer independence 
(no flux of particles or gases entering or leaving an individual layer), special 
provision must be made for spatial changes in the horizontal meteorology and 
for gravitational settling or precipitation scavenging. In addition, the type of 
source within a layer must be considered; that is, whether there is a ground 
cloud source or a plume cloud source. 

The NASA/MSFC Multilayer Diffusion Model has six models (Fig. 2) 
which account for three categories of dispersive constraints: the source dis- 
tribution, the environmental effects, and the depositional effects. This flexi- 
bility is required to deal with the stages of the development of the exhaust cloud 
and the complex, potentially varying, meteorological conditions. These models 
can be used alone to describe all the environmental layers or in superimposed 
combinations where variations in layer meteorology require different modeling. 
For the Introductory overview, however, these combinations will not be con- 
sidered. The primary output of all models is a mapping of the regimes of the 
concentration and dosage isopleths and centerline profiles for concentration 
and dosage. 
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The fundamental category of dispersive constraints is the source distri- 
bution. The two distributions ares 

1. The elliptic-cylindrical source which assumes a two-dimensional 
Gaussian distribution in the x-y plane and a uniform distribution in the vertical 
direction. 


2. The ellipsoidal source which assumes a three-dimensional Gaussian 
distribution. 

Model 1 is for the elliptic-cylindrical source whose vertical expansion 
is constrained by the layer boundaries - thus it has only a two-dimensional 
expansion in the horizontal plane due to turbulence mising. This model is 
normally used to describe the rocket’ s inflight plume cloud. 

Model 3 is for the ellipsoidal source and is assumed to expand in all 
three dimensions as the effluents are propagated downstream. When the 
ellipsoidal source reaches the top of the mixing layer, the distribution of the 
constituents is reflected back into the expanding vertical distribution. On the 
other hand, that fraction not lost in surface deposition is also reflected back 
in a similar manner. After sufficient mixing, the ellipsoidal distribution 
becomes an elliptic-cylindrical distribution (Model l). While Model 3 is nor- 
mally used to describe the dispersion of the rocket’ s ground cloud, it could be 
used to model upper air explosions. The formulation for Model 3 has been 
provided in Appendix B. 

The second category is environmental effects. The two effects are: 

1. No turbulence mixing in the upper atmosphere. 

2. Changes in meteorological conditions as the constituents are trans- 
ported downstream. 

Model 2 is the same as Model 1 except it is assumed that there is no 
turbulent mixing. This implies that the exhaust material Just meanders along 
the layer without dispersing. While Model 2 is not r nerally used, movies 
of rocket firings clearly show that under some special meteorological conditions 
this model is required. While the Multilayer Diffusion Model is general in 
applicability, it is specific in meteorological parameters and launch description. 

Model 4 updates the diffusion model with changes in meteorological con- 
ditions and structure which can occur as the constituents propagate downstream. 
This model assumes that the vertical concentration of material has become 
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uniform throughout each layer when a step-change in the meteorological con- 
ditions is introduced, resulting in the destruction of the original layer bound- 
aries and the formation of new layer boundaries. The concentration fields 
which exist at this time are treated as new sources. In those new layers which 
now comprise more than one old layer, the old concentration is mapped as two 
independent concentration sources and then superimposed for the resulting con- 
centration and dosage mappings. 

The third category of dispersive contraints includes the deposition due 
to: (l) precipitation scavenging and (2) gravitational settling. 

Model 5 accounts for precipitation scavenging. An example of where 
Model 5 must be used is in solid rocket launches during the occurrence of rain, 
because the HC1 will be scavenged by the rain. Model 6 describes the ground 
deposition due to gravitational settling of particles or droplets. Wind shears 
are incorporated in this model to account for the effect of the settling velocity 
of the particulate matter. There are two forms for the source in this model; 
namely, 

1. The source that extends vertically through the entire layer with a 
uniform distribution — this is the same source model as used with Models 1 
and 2. 


2. A volume in the Kth layer — this is the same source model as used 
with Model 3. 

Model 6 is very important in the analysis of the settling of AljOj particles 
released in solid rocket firings. 

The treatment of cold spills and fuel leaks that occur near ground level 
requires a continuous source, but the models that have been considered so far 
are for discrete sources; therefore, the models must be adapted for the use in 
predicting concentration-dosage levels downwind from continuous sources. 

The layer of the environment influenced by the ground-level spills and 
leaks can be treated as homogeneous; therefore, the general formulas for con- 
centration and dosage f equations (14) and (15)] presented in the initial discus- 
sion would be applicable if spills and leaks are treated as continuous sources. 

To visualize this adaption for these formulas, assume a source cloud with a 
concentration distribution that implies a gtven dosage at a point for this cloud; 
that is, the dosage per event. If there are a number of similar clouds, dis- 
cretely spaced, then for each cloud we obtain a dosage for each cloud whose 
sum corresponds to the total dosage for the entire event. 
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In summary, the Multilayer Diffusion Model is composed of six sc 1 '- 
models. Models 1 and 3 are designed to distinguish between the two sources 
of toxL' cloud formation — the ground cloud during the initial launch phase 
(Model 3) and the plume cloud after the Initial launch phase (Model l). From 
the standpoint of environmental impact, the description of the fields of the 
ground deposition of materials from the ground cloud is of primary significance; 
this description is afforded by Model 3. Generally, this model is employed in 
the surface layer, but can be employed in any layer where the source ooes not 
extend through the entire layer. 

Model 2 was designed to account for a lack of turbulent mixing which can 
occur in the upper atmosphere. Model 4 is employed when a change in meteor- 
ological condition occurs during the downstream transport of the cloud. In the 
event of rain, the precipitation scavenging of both gases and particles can be 
accounted for in Model 5. The fallout of particulate matter on the ground is 
the domain of Model 6. These six submodels form the basic set of equations 
which are available to treat the diffusion problem. To model a specific launch 
of a vehicle, it is necessary to blend these equations and adjust the model 
parameters to the specific meteorological conditions of the launch, to the 
specific terrain around the launch site, and to the specific vehicle being 
launched; thus the degree of complexity in the diffusion model. 


RESULTS OF THE CALCULATIONS 


The results of the model concentration computations and descriptions of 
the meteorological regimes at Kennedy Space Center and Vandenberg Air Force 
Base, for which the calculations were maae, are as follows. 

Peak ground-level and 10-min average concentration calculations were 
made for the nine meteorological regimes at KSC and the six meteorological 
regimes at VAFB as indicated in Table 2. As noted in the table, the first three 
regimes listed under KSC are "typical” regimes. Also in Table 2, ground-level 
concentrations were calculated for a meteorological regime in which a stationary 
front extended from east to west through central Florida in the vicinity of KSC 
on October 2, 1972. Ground-level concentrations were also calculated for three 
meteorological regimes associa+^d with the approach and passage of a cold front 
at KSC during the period from October 19 to 21, 1972. Calculations were made 
for another case in which a cold front was located south of KSC on November 26, 
1972. Finally, concentrations on November 27, 1972, were calculated for a 
regime where a fair-weather, high pressure regime existed over the South- 
eastern United States. 
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As presented In Table 2, concentration calculations were made at VAFB 
for three meteorological regimes in which the frontal activity influences dis- 
persion. Calculations were made for October 10, 1972, (stationery upper- 
level trough west of VAFB) , January 16, 1973, (cold front northwest of VAFB), 
and January 17, 1973, (cold front south of VAFB). 


Concentration Calculations for Kennedy Space t :?r 

Fall, Spring, and Sea-Breeze Meteorological Regimes . The meteoro- 
logical data used in the concentration calculations for the fall, spring, and sea- 
breeze regimes at KSC were derived from the mean monthly wind speed, wind 
direction, and temperature profiles for KSC [ 22, 23] . These profiles have 
been used in previous hazard calculations for launches at KSC f 23] . A study 
of the KSC climatology indicated that the average depth of the surface mixing 
layer in the fall season associated with the easterly winds required to transport 
the ground cloud inland is approximately 1000 m. During the spring, there are 
occasions when the surface mixing layer reaches a depth of 2000 m. The after- 
noon sea breeze, which Is common to all seasons, has an average surface mix- 
ing layer of 300 m. The vertical profiles of wind speed, wind direction, and 
temperature for these three regimes are shown in Figures 3, 4, and 5. Table 

3 gives the calculated HC1, CO, and A1 2 0 3 peak concentrations downwind from 
the normal launch of a Space Shuttle vehicle for the three regimes at KSC. The 
slow decrease of ground-level concentrations with distance is evident for all 
three regimes at some distance between 1 and 10 km and reflects the gradual 
change in the vertical distribution of material from the initial Gaussian form 

to a rectangular form brought about by multiple reflection between the base and 
top of the surface layer. The highest concentration in these three cases was 
the fall regime. Outside the KSC industrial area of 10 km, HC1 was less than 
2. 0 ppm, 3. 5 for CO, and 4. 2 mg/m s for A1 2 0 3 as presented in Table 3. Table 

4 gives the 10-min average concentration for HC1 at KSC ranging from 1. 39 to 
less than 0. 01 ppm from 1 to 80 km for all nine meteorological conditions 
listed, which is well below the applicable limits developed by the NAS/NRC 
Committee on Toxicology f 24, 25] . 

Pre-Cold Front, Cold Front, and Post-Cold Front Meteorological 
Regimes . Ground-level concentrations downwind from the launch pad were 
calculated for three meteorological regimes associated with the approach and 
passage of a cold front. The meteorological data of October 19, 20, and 21, 
1972, vere selected to be representative of these regimes at KSC. Figures 6, 
7, and 8 show the vertical profiles of wind speed, wind dir action, and tempera- 
ture at 0700 Eastern Standard Time (EST) for the three days. At 0700 EST on 
October 19, a cold front was located northwest of Florida, extending from 
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Georgia through Southeast Alabama and then westward. Florid x weather was 
fair with scattered high clouds and local haze conditions. The KSC 0700 EST 
sounding showed a stable layer from the surface to about 21 1 m above the 
ground. At 0700 EST on October 20, the cold front was oriented from east to 
west and was located just to the south of KSC. Rain and rain showers were 
occurring in the vicinity of the front. The KSC 0700 EST sounding Indicates 
the presence of a moist, unstable air mass over KSC with a r irface mixing 
layer extending to about 2000 m. The cold front had moved south into the Straits 
of Florida by 0700 EST on October 21. The 070C EST sounding at KSC indicates 
a surface mixing layer depth of approximately 1400 m dth dry, warm air aloft. 

The calculated HC1 peak ground-level concentrations as given in Table 3 
at a downwind distance of 10 km are 0. 7, 1. 0, and 2. 6 ppm for the pre-cold 
front, cold front, and post-cold front at KSC. 

Additional Meteorological Regimes at KSC . Table 3 gives the peak con- 
centration calculations for a stationary front south of KSC, a cold front south 
of KSC, and a fair-weather, high pressure regime at KSC for measurements at 
10 km on October 2, 1972, November 26, 1972, and November 27, 1972, 

Peak concentration calculations for HC1 at the 10-km downwind distance 
were 1.2, 3.2, and 0-2 ppm, respectively, as presented in Table 3. Table 4 
gives the 10-min average concentration for HC1 at KSC ranging from 1. 39 to 
0. 01 ipm or less for all the meteorological regimes from 1 to 80 km, which 
is well below the applicable exposure limits of 4 ppm developed by the NAS/NRC 
Committee on Toxicology. 


Concentration Calculations for Vandenberg Air Force Base 

Morning and Afternoon Fair-Weather Regimes, hi the absence of frontal 
activity, the weather in the lower 2000 to 3000 m at VAFB is dominated by the 
land-sea-breeze regime. The marine (subsidence) inversion is present at 
VAFB over 90 percent of the time in the summer and approximately 50 percent 
of the time in the winter. With the Pacific high centered to the west and a 
thermal low common over the California interior, the prevailing gradient wind 
in all seasons is from the north or notheast. Cold air drainage from canyons 
in the vicinity of the launch areas contribute to the offshore winds at night and 
during the morning. Figure 9 shows typical morning profiles of wind speed, 
wind direction, and temperature during this type of fair-weather regime. At 
the time of the sounding, the westerly sea breeze is beginning to develop near 
the surface with northeasterly winds above. The surface mixing layer extends 
to the base of the inversion which is about 400 m above the ground. 
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During the day the layer of onshore flow deepens as the sea breeze 
becomes well established. Figure 10 shows typical wind speed, wind direction, 
and temperature profiles for an afternoon sea-breeze regime; the inversion 
base is about 225 m above the surface. Relatively high inversions are also 
observed at VAFB, often associated with an upper level trough to the west. 
Figure 11 shows typical wind and temperature profiles for an afternoon sea- 
breeze regirr with an inversion base approximately 775 m above the surface. 

The results of the peak concentration calculations for the morning low 
level meteorological regime at VAFB are presented in Table 3. In Table 3 the 
peak concentration profiles for the afternoon sea-breeze, low level meteoro- 
logical regime are given. At 10 km, the morning inversion measured 1.70 
ppm, and the afternoon inversion was 2.4 ppm. Also as given in Table 3 the 
sea-breeze, high level inversion meteorological regime measured approxi- 
mately 6. 5 ppm at 10 km. 

Additional M eteorological Regimes at VAFB. Table 3 presents the peak 
concentration calculations for a stationary, upper level trough west of VAFB, a 
cold front northwest c' VAFB, and a cold front south of VAFB. HC1 concentra- 
tio " at 10 km were about 1. 5, 1. 2, and 1. 8 ppm. 

Table 5 shows the 10-min average concentrations for HC1 at VAFB, 
ranging from 1. 75 to 0. 01 ppm from 1 to 80 km for all the meteorological 
regimes. These values are well below the applicable exposure limits of 4 
ppm developed by the NAS/NRC Committee on Toxicology. 


Summary of Results 

Table 3 lists the calculated peak ground-level concentrations for each 
of the three pollutants of interest at downwind distances of 1, 2, 5, and 10 km 
from the launch pad for normal launches of the Space Shuttle at Kennedy Space 
Center and Vandenberg Air Force Base. The pollutant concentrations at dis- 
tances of 5 and 10 km are enclosed by solid lines to indicate that they should 
receive primary attention in assessing the potential environment hazard of the 
Space Shuttle emissions. The model concentrations made concern the vertical 
distributions of material in the stabilized cloud of exhaust products. Accurate 
measurements of this distribution ar ot available. Evidence from photo- 
graphic data and very limited ground-level pollutant concentration measure- 
ments strongly indicates that only very small concentrations of exhaust products 
are found near ground-level after the cloud has stabilized. High concentrations 
of exhaust products exist in the immediate vicinity of the launch pad during the 
first minute or so after ignition as the result of the lateral deflection of hot 
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exhaust products which may extend outward from the pad to a maximum radial 
distance of about 1 km. However, most of the hot exhaust products then ascend 
and are incorporated into the main body of the stabilized ground cloud. In our 
model calculations the assumption is that the exhaust products in the stabilized 
cloud have a Gaussian distribution which is centered at the cloud stabilization 
height and extends to the ground. This assumption appears to result in larger 
ground-level concentrations near the launch pad after cloud stabilization than 
actually occur, which in turn leads to overestimates of downwind ground-level 
concentration s in the first few kilometers. 

The calculated peak ground-level concentrations at a downwind distance 
of 10 km are appropriate for use in assessing the potential hazard to uncon- 
trolled populations because, at both KSC and VAFB, this is the minimum dis- 
tance from the launch pad area to the boundaries of uncontrolled population 
areas. Under this assumption, the calculated HC1 peak ground-level concentra- 
tions in Table 3 at KSC for normal Space Shuttle launches range from 3. 2 to 
0. 2 ppm, and the HC1 peak ground-level concentrations at VAFB range from 
6. 5 to 1. 2 ppm. 

As shown in Tables 4 and 5, the 10-min average concentrations for HC1 
lange from 1. 39 to less than 0. 01 ppm from 1 to 80 km at KSC. At VAFB 
these values range from 1. 75 to 0. 01 ppm which is well below the applicable 
exposure limits of 4 ppm developed by the NAS/NRC Committee on Toxicology. 

The 10-min average concentrations for CO for all the meteorological 
conditions (Tables 4 and 5) at KSC were 2.45 to 0.01 ppm from downwind dis- 
tances 1 to 80 km from the launch pad, and at VAFB the measurements were 
3. 08 ppm to 0. 02 ppm. These measurements for Aip 3 were 3.03 mg/m 3 to 
0. 01 mg/ m 3 for KSC, and 3. 97 to 0. 03 mg/m 3 for the same conditions at VAFB. 


18 



TABLE 1. FUEL PROPERTIES OF SPACE SHUTTLE ENGINES 


Total Mass Flow Rate (g/sec) 


All Solid Engines 

<w s ) 

9.45126 x 10 6 

Liquid Engines 

(W L ) 

1. 53316 x 10 6 

Fuel Heat Content (cal/g) 



Solid Engines 

<H s> 

691 

Liquid Engines 

<v 

500 

Fuel Composition (percent by weight) 


HC1 

(FM) 

20.7 

CO 

(FM) 

28.0 

AljOj 

(FM) 

30.4 




TABLE 2. METEOROLOGICAL REGIMES FOR WHICH CONCENTRATION 

CALCULATIONS WERE MADE 


Location 

Meteorological Regime 

Date 

KSC 

Fall 

Spring 

Sea-Breeze 

Typical Regimes 


Cold Front North of KSC 

October 19, 1972 


Cold Front Near KSC 

October 20, 1972 


Cold Front South of KSC 

October 21, 1972 


Stationary Front South of KSC 

October 2, 1972 


Cold Front South of KSC 

November 26, 1972 


Fair Weather, High Pressure 

November 27, 1972 

VAFB 

Morning 

Sea-Breeze with Low Level Inversion 
Sea-Breeze with High Level Inversion 

Typical Regimes 


Stationary Upper-Level Trough West of VAFB 

October 10, 1972 


Cold Front Northwest of VAFB 

January 16, 1 973 


Cold Front South of VAFB 

January 17, 1973 
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TABLE 3. SUMMARY OF PEAK GROUND-LEVEL CONCENTRATIONS 
FOR NORMAL LAUNCHES OF THE SPACE SHUTTLE AT KSC AND 
VAFB UNDER VARIOUS METEOROLOGICAL CONDITIONS 


Meteorological 

Regime 

Pollutant 

Centerline Ground-Level 
Concentration 

Distance from Launch Pad (km) 

1 

2 

5 

lOl 

a. Kennedy Space Center 

Fall 

HC1 (ppm) 

mm 


4.4 

2.0 


CO (ppm) 


8.0 

7.5 

3.5 


A1 2 0 3 (mg/m 3 ) 


10.1 

9.2 

4.2 

Spring 

HC1 (ppm) 

4.6 

2.4 

2.0 

1.6 


CO (ppm) 

8.2 

3.6 

3.2 

3.0 


Al 2 °3 (mg/m 3 ) 

10.0 

5.2 

4.4 

3.6 

Sea Breeze 

HC1 (ppm) 

14.0 

7.0 

1.6 

■a 


CO (ppm) 

24.0 

12.5 

2.8 



AljOj ( mg/m 3 ) 

32.2 

15.5 


■b 

Cold Front North of 

HC1 (ppm) 

19.0 

7.0 


0.7 

KSC (Oct. 19, 1972) 

CO (ppm) 

33.0 

12.5 

4.0 

1.3 


A1jQ 3 (mg/m 3 ) 

42.0 

16.0 


■a 

Cold Front Near 

HC1 (ppm) 

50.1 

20.0 

1 

1.0 

KSC (Oct. 20, 1972) 

CO (ppm) 

90.0 

34.0 

1 



Aip 3 (mg/m 3 ) 

100.0 

44.0 



Cold Front South of 

HC1 (ppm) 

9.2 

3.6 



KSC (Oct. 21, 1972) 

CO (ppm) 

15.0 


5.0 



Aip 3 (mg/m 3 ) 

24.0 

8.0 

7.0 

6.0 

Stationary Front 

HC1 (ppm) 

3.2 



1.2 

South of KSC 

CO (ppm) 

5.5 



2.1 

(Oct. 2, 1972) 

A1j 0 3 (mg/m 3 ) 

7.0 

■m 

5.0 

2.7 
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TABLE 3. (Concluded) 


Meteorological 

Regime 

Pollutant 

Centerline Ground-Level 
Concentration 

Distance from Launch Pad (km) 


2 

5 

10 

Cold Front South of 

HC1 (ppm) 

10.1 

4.2 

3.5 

3.2 

KSC (Nov. 26, 1972) 

CO (ppm) 

20.0 

7.5 

6.0 

5. 8 


Aip 3 (mg/m 3 ) 

30.0 

9.8 

7.5 

7.5 

Fair Weather, High 

HC1 (ppm) 

9.5 

4.4 

0. 8 

0.2 

Pressure 

CO (ppm) 

16.0 

7.5 

1.5 

0.4 

(Nov. 27, 1972) 

Aip 3 (mg/m 3 ) 

20.0 

10.0 

2.0 

0.5 

b. Vandenberg Air Force Base 

Morning Inversion 

HC1 (ppm) 

20.0 

11.5 

5.5 

1.7 


CO (ppm) 

34.0 

19.7 

9.5 

3.0 


Aip 3 (mg/m 3 ) 

45.0 

24.0 

12.0 

3.8 

Sea Breeze with Low 

HC1 (ppm) 

,21. 1 

11.2 

4.7 

2.4 

Level Inversion 

CO (ppm) 

37.2 

19.7 

8. 2 

4.2 


A1 2 0 3 (mg/m 3 ) 

47.6 

25.3 

10.5 

5.3 

Sea Breeze with 

HC1 (ppm) 

17.3 

11.7 

10. 5 

6.5 

High Level Inversion 

CO (ppm) 

30.4 

20.5 

18.5 

11.4 


Aip 3 (mg/m 3 ) 

39.2 

26.5 

23. 8 

14.7 

Stationary Upper 

HC1 (ppm) 

5.4 

2.3 

1.9 

1.5 

Level Trough West of 

CO (ppm) 

9.5 

4.1 

3.5 

2.7 

VAFB (Oct. 10, 1972) 

Aip 3 (mg/m 3 ) 

12.0 

5.4 

4.4 

3.5 

Cold Front Northwest 

HC1 (ppm) 

17.0 

7.0 

3.6 

1.2 

of VAFB 

CO (ppm) 

28.0 

12.0 

6.2 

2.0 

(Jan. 16, 1973) 

Aip 3 (mg/m 3 ) 

36.0 

16.0 

8.0 

2.7 

Cold Front South 

HC1 (ppm) 

7.0 

2.6 

1. 8 

1.8 

of VAFB 

CO (ppm) 

8.2 

3.2 

3.0 

3.2 

(Jan. 17, 1973) 

Aip 3 (mg/m 3 ) 

16.0 

i 

6.0 

4.0 

4.0 




TABLE 4. SUMMARY OF 10-min AVERAGE CENTERLINE 
CONCENTRATIONS FOR SPACE SHUTTLE LAUNCHES AT 
KSC (POLLUTANT CONCENTRATIONS FOR HC1 AND CO 
ARE IN ppm WHILE THOSE FOR A1 2 0 3 ARE IN mg/m 3 ) 


Meteorological 

Regime 

— 

Pollutant 

Centerline Ground- Level 
Concentration 

bistance from Launch Pad (km) 

1 

2 

5 

10 

20 

40 

80 

Fall 

HCl 


0.38 

0.69 

M 

0.22 

0.18 



CO 

1 3 

0.66 

1,22 

Egg 

0.38 

0.19 



A1 2 0 3 

0.54 

0. 82 

1, 51 

m 

0.48 

0. 23 

0.09 

Spring 

HCl 

0.10 

0.10 

0.25 

0. 29 

0.18 

0.09 



CO 

0.16 

0.14 

0.42 

0.52 

0.33 

0.16 



A1 2 0 3 

0.22 

0.21 

0.54 

0.64 


0.20 

0.15 

Sea Breeze 

HCl 

0.42 

0.38 

0.19 

0.10 

0.05 

0.02 

0.01 


CO 

0. 73 

0.67 

0.33 

0.17 


0.04 

0.01 


Alp3 

0.92 

0.84 

0.42 

0.22 

0.11 

0.05 

0.02 

Stationary Front 

HCl 

0.38 

0.60 

,1.39i 


0.33 

0.09 

0.02 

South of KSC 

CO 

0.66 

1.06 

2.45 

1,69 

0.59 

0.16 

0.04 

(10-2-72) 

Aip 3 

0.82 

1.72 

3. 03 

2.10 

0.73 


0.05 

Cold Front North 

HCl 

0.39 

0.30 

0.20 

0.10 

na 

0. 02 

0.01 

of KSC 

CO 

0.68 

0. 53 

0.34 

0.18 


0.04 

0. 02 

(10-19-72) 

Aip 3 

0.87 

0.68 

0.44 

0.23 

0.11 

0. 06 

0.02 

Cold Front Near 

HCl 

0.57 

0.37 

0.17 

0.09 


0.02 

0.01 

KSC (10-20-72) 

CO 

1.01 

0.64 

0.31 

0.16 

0.08 

0.04 

0.02 


Alp 3 

1.26 

0. 81 

0.38 

0. 20 

0.10 

0.05 

0.03 

Cold Front 

HCl 

0.17 

0.15 

0.37 

0.51 

0.32 

0.16 

0.06 

South of KSC 

CO 

0.88 

0.87 

0.80 

0.94 

0. 57 

0.29 

0.11 

(10-21-72) 

Aip 3 

0.37 

0.33 

0.81 

1.14 

0.71 

0. 36 

0.14 

Cold Front 

HCl 

0.25 

0.19 

0.39 

0.71 


0.28 

0.09 

South of KSC 

CO 

0.43 

0.33 

0.68: 

lr 25 

Bps 

0.49 

0.16 

(11-26-72) 

A1 2 0 3 

| 0.55 

0.42 

0.86 

1.59 

1.34 

0,63 

0.20 

Fair Weather 

HCl 

0.62 

0.47 

0.19 

0.10 

0.04 

© 

0 

0,00 

High Pressure 

CO 

1.09 

0.83 

0.34 

0.17 

0.07 

0.02 

0.01 

(11-27-72) 

Aip 3 

1.46 

1.10 

0.46 

0.22 

0.10 

0.03 

0.01 
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TABLE 5. SUMMARY OF 10-min AVERAGE CENTERLINE CON- 
CENTRATIONS FOR SPACE SHUTTLE LAUNCHES AT VAFB 
( POLLUTANT CONCENTRATIONS FOR HC1 AND CO ARE IN 
ppm WHILE THOSE FOR A1 2 0 3 ARE IN mg/m 3 ) 


Meteorological 

Regime 

Pollutan* 

C enterl ine Ground- Level 
Concentration 

Distance from Launch Pad (km) 

1 

2 

5 ' 

10 

20 

40 

80 

Morning 

HC1 

0.69 

0. 81 

0.79 

0.42 

0.21 

0.10 

0.03 


CO 

1.10 

1.30 

1.34 

0.72 

0.37 

0.15 

0.05 


Aljp 3 

1.43 

1.69 

1.74 

0.94 

0.48 

0.20 

0.06 

Sea Breeze with 

HC1 

0.63 

0.46 

0.20 

0.10 

0.05 

0.07 

0.01 

Low Level Inversion 

CO 

1.11 

0. 81 

0.35 

0.18 

0.08 

0.05 

0.02 


AIjPs 

1.42 

1.04 

0.45 

0. 23 

0.12 

0.07 

0.03 

Sea Breeze with 

HC1 

0.63 

0.85 

1.75 

1.18 

0.61 

0.31 

0.16 

High Level Inversion 

CO 

1.11 

1.50 

3.08 

2. 07 

1.08 

0.55 

0.28 


AljOj 

1.43 

1.93 

3.97 

2.67 

1.39 

0.71 

0.36 

Stationary Upper 

HC1 
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Figure 2. Block diagram of the computer program for the 
NASA/MSFC Multilayer Diffusion Model. 
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Figure 5. Vertical profiles of temperature, wind speed, and wind direction 
for the sea-breeze meteorological regime atKSC. 



























Figure 11. Vertical profiles of temperature, wind speed, and wind direction for the afternoon 
sea-breeze, high level inversion meteorological regime at VAFB. 






APPENDIX A 


CLOUD RISE FORMULAE [l] 


In order to determine whether an adiabatic or stable cloud rise relation 
should be utilized, it is necessary to determine the vertical potential tempera- 
ture gradient (90/9z) which is described by 


90 e 9T e 9T g /A , 

— = — + - 6 - « + - 6 - ( A—] 

9z T 9z c 9z c y 

P P 

and g/c is equal to the adiabatic lapse rate (10° C/km). 

P 

In addition, 0 and T are the potential and ambient temperature, g is 
the gravitational acceleration, and c is the specific heat of air. If 


r * o 

9z 


the adiabatic cloud rise relation is used. 


(A-2) 


The maximum cloud rise z downwind from an instantaneous source 

ml 

in adiabatic atmosphere is given by 



(A-3) 


whereas, the maximum cloud rise 
in a stable atmosphere is given by 


z , downwind from an instantaneous source 
ml 



(A-4) 
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where is the instantaneous buoyancy parameter (3gQj/47rpCpT) , is 
the effective heat released, p is the density of ambient air, y^ is the entrain- 
ment coefficient, r is the initial cloud radius at the surface, s accounts for 
K 

the vertical gradient of the potential temperature, and t g j is the time to reach 

stabilization. The subscript I means instantaneous and is used to flag a dif- 
ference in the cloud rise models. The buoyancy terms, which are a function of 
the heat released and the type of entrainment, spherical and cylindrical, reflect 
the major difference in the two sources. 

Equations (A-3) and (A-4) assume that the initial upward momentum 
imparted to the exhaust gases by reflection from the ground surface and launch 
pad hardware is insignificant in comparison with the diermal buoyancy flux. 
These relations are normally used with solid rocket motors. 


The following formulae for the maximum buoyant rise of clouds from 
continuous sources are also based on procedures similar to those given by 
Briggs ( 1970). The maximum cloud rise downwind from a continuous 

source in an adiabatic atmosphere is given by 


z 

me 


3F x 2 
c sc 


2y 2 u s 
c 


+ 




(A-5) 


The maximum cloud rise z 
stable atmosphere is given by 


downwind from a continuous source in a 



(A-6) 


where F is the continuous buoyancy flux parameter and is equal to 
c 

gQ /it pc T. The subscript c implies that the associated parameter is unique 
c p 

to the continuous source. The primary difference in these continuous source 
relations is that the temperature constraint in the stable atmosphere results 
in a buoyancy damping. 
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Equations (A-5) and (A-6) assume that the initial momentum flux 
inparted to the cloud by dynamic forces is negligible in comparison to buoyancy 
flux. Again, experience in calculating cloud rise for normal launches of large 
liquid fueled rockets and for static firings has shown that this assumption is 
reasonable. 
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APPENDIX B 


CONCENTRATION-DOSAGE FORMULATION FOR NASA/MSFC 
MULTI LAYER DIFFUSION MODEL 


The fundamental relation for the concentration-dosage calculation will 
be presented for the ellipsoidal source used in Model 3. These relations are 
appropriate to the elliptic-cylindrical distribution of Model 1 if the vertical 
dispersive interaction is neglected. This part of the appendix is complex; 
therefore, its use is recommended only when a detailed scientific knowledge 
is required. 

The dosage equation for Model 3 in the Kth layer is given by the 
expression 


°K { x K’ y K> z BK <Z K < Z TK } = 


Q K I ' y K. a 

— ( exp 

2irff yK a zK U K I 2a yK 


exp 


< H K- z K> a 


2 "zK 


exp I 


( H K - 2z 3K + Z K> 2 
2a zR 2 


+ y | exp | ' ^ 2 ‘ (Ztk ' Zbk) ' (Hr ~ 2zbk * Zr) ) 2 


R 2 * 

I" ( 2 '( z TK ■ Z BK) + («K • z k)) 1 1 

L 2 ^ J 

| ( 2i ( z TIC ' Z BK> ~ (»K ' Z K)) 2 1 

+ exP L 2o * J 


i=l 


+ exp 


+ exp J '^ 2 ‘ (ZTK ' ZBK) 1 (HfC ‘ 2ZbK 1 ZK ^ 2 

~ a zK 


2 a ,-- 1 


(B-l) 
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where Q„ corresponds to the source strength or total mass of material in the 
K 

layer and H is the height of the centroid of the stabilized cloud. 

The standard deviation of the vertical dosage distribution (a^) is 
defined by the expression 


zK 


= a* 


*K 


+ x 


zK 


X rzK^ 


EK rzK 


X rzK 


(B-2) 


where < 7 * describes the mean standard deviation of the wind elevation angle, 

gives the vertical virtual distance, accounts for vertical diffusion, 

and x _ is the distance over which rectilinear vertical expansion occurs 
rzK 

downwind from an ideal point source in the Kth layer. 


In the surface layer (K = l), the standard deviation of the wind ele- 
vation angle (a RR ) at the height z R is described by 


EK 


(k = 1} = 


ER 


(z tk (K= 1 }) 


(q+ D(z tk (K= 1} -z R )(z R ) 


»vq+l / ^+1 

L—li!?! 1 (JL) 

\ 18 °/ ’ 


(B-3) 


where the power-law exponent (q) for the vertical profile of the standard 
deviation of the wind elevation angle in the surface layer is 



(B-4) 


where a { K = l} is the standard deviation of the wind elevation angle at 
ETK 

the top of the surface layer. Above the surface layer (K > l), the standard 
deviation of the wind elevation angle is 
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4 


EK {K>1) ^ETf; + a EBK^ (360) 


(B-5) 


where c r , 


and cr . 


ETK EBK 

angle at the top and the base of the layer. 


are the standard deviations of the wind elevation 


The vertical virtual distance x is given by the expression 


c ( K ) 

zo 


- X 


EK 


RzK 




X rzK la* 


zo 


EK X rzK 


<j {K} s o' x 
zp EK rzK 


Up 


K 


- x + x V (1 - (3 ) ; <7 {K} 2= <r* x 

RzK rzK K zo 


EK rzK 

(B-6) 


where a { K} is the standard deviation of the vertical dosage distribution at 
zo 

x„ _ _ , the distance from the source where the measurement is made in the Kth 
RzK 

layer. 

The remaining terms are common also to Model 1; that is, what has just 
been discussed is to account for the vertical expansion of the source cloud. 

The quantity u^ in equation (B-l) is the mean cloud transport speed 

in the Kth layer. In the surface layer (K = l) , the wind speed-height profile 
is defined according to the power-law expression 


u(z K * K= l) = u R 


z k (K = 1} 


(B-7) 


where u_ is the mean wind speed measured at the reference height z_ and 

K x\ 

the power-law exponent (p) for the wind speed profile In the surface layer is 
described by 
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p = 


log 


u tk^ k " ^ 


u 


R 





(B-8) 


Here, u^{ K = l) corresponds to the mean wind speed at the top of the surface 
layer (z^{K = l) ). Thus, In the surface layer, the mean cloud transport 
speed (u(K = l) ) is 




“r 

< Z TK {K = * Z R> Z R P 


TK 


f (z K {K= l}) P dz , (B-9) 


which reduces to 


yK= 1} = 


u 


R 


< Z ™{K-») ltP -(V 1+P 


(B-10) 


(z im { K = 1} - z R )(z R ) P (1+ p) 


In layers above the surface layer (K > l) , the wind speed-height profile 
(u{ z R , K > l) ) is assumed linear and is defined as 


G{z k ,k>i) = G bk + 


U TK ~ U BK \ 
Z TK " Z BK / 



z 



(B-ll) 


where and u RR describe the mean wind speed at the top of the layer and 

at the base of the layer, respectively. In the Kth layer (K >1) , the mean 
cloud transport speed (u^{K>l|) is 


y K > 1} - (J ra +=B K )/2 


(B-12) 
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The standard deviation of the crosswind dosage distribution 
defined by 


(< V 


is 


yK 


a 


'akW x ryK i V> . _ ~ V 


K 


O',, x 
K ryK 


M K*1{ 


4.3 


*/* 


(B-13) 


where cr^{ t^} corresponds to the mean layer standard deviation of the wind 

azimuth for the cloud stabilization time (r ). In the surface layer (K=l), 

K 


. , „ L'“TK 

< 'ak { V k = i) = 


( I „{K = l}) ra * 1 -(z R )^ 1 


( m+ 1 )(z tk (k= 1} - z R )(z R ) 


m 


(B-14) 


where the standard deviation of the wind azimuth angle (o* {t r }) at height 

z and for the cloud stabilization time t is 
H K 


r * { T 

AR l K 


.} « 


r , / T K\ 1/5 / 7T \ 

' AR oK> V t „k/ V 180 / 


(B-15) 


Here <7 An { r } is the standard deviation of the wind azimuth angle at height 
Ait oK. 

z R and for the reference time period (t ) , and the power-law exponent (m) 

for the vertical profile of the standard deviation of the wind azimuth angle in 
the surface layer is 



(B-16) 



Then, 



where » ATK {V’ K ’ 1 > 


is the standard deviation of the wind azimuth angle 


at the top of the surface layer for the reference time period. For layers above 
the surface (K > l), 


^ATK^'K* K > ^ 


^ATK^K^ 



where, 


ATK K 


( T J = * J 


ATK oK’ 



(B-18) 


(B-19) 


Here a . r } is the standard deviation of the wind azimuth angle at the 
ATK oK 

top of the layer. 


r» { T } 

ABK X K f 


r { T 1 

ABK X oK' 



(B-20) 


where 7 {r } is the standard deviation of the wind azimuth angle in 
ABK OK 

degrees at the base of the layer for the reference time period (t ._) , 

OK 
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< — 


The crosswind virtual distance is 


<r S K) 


> ■ ■ **»* 


(B-21) 


when 


°yo iK) S °"ak {t K } x ryK 


a K X ryK 


<7 {K} 

VO 

< ’ak {t K> X ryK 


when 


- x RyK + x ryK (1 ' V 


(B-22) 


a yo^ " a AK^ T K^ X ryK 


Here, a {k} is the standard deviation of the lateral source dimension in the 
yo 

laver at downwind distance x , x is the distance over which rectilinear 

KyK ryK 

crosswind expansion occurs downwind from an ideal point source, and a 

K 

describes the lateral diffusion in the layer. The vertical wind direction shear 

(A O') in the layer is 
K 


A0 k ^TK -<9 BK^ (l8o) 


(B-23) 


where 9 and 6 are the mean wind direction at the top and at the base 
IK BK 

of the layer, respectively. 


I 


'L.__ - 


The concentration algorithm is of the same form for the first three 
models; however, the dosage term (D ) does depend on which model has 

been utilized, and thus adjusts the concentration description to the specific 
model of interest. 

The maximum concentration for the first three models in the Kth 
layer is given by the expression 


X K { V y K> Z K } * 


d k"k 


(B-24) 


\/~2ir <t 


xK 


where the standard deviation of the alongwind concentration distribution (ct^) 
in the layer is 


xK 



1 V 2 


4.3 


+ <7 2 {K} 

xo 


(B-25) 


and the alongwind cloud length (L{x^}) for a point source in the layer at the 
distance x^ from the source is 


l{x k } ■ 


0.28 (A5 k )(x r ) 


0 


AUr^O 


(B-26) 


where Au^, is the vertical wind speed shear in the layer and is defined as 


A^{K =1} = u tk {K= 1} - u. 


R 


(B-27) 
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or 


Ay K >n = 



(B-28) 


and a { K} is the standard deviation of the alongwind source dimension in 
xo 

the layer at the point of cloud stabilization. The above equation for L{ x^} 

is based on the theoretical and empirical results reported by Tyldesley and 
Wallington f 26] who analyzed ground-level concentration measurements made 
at a distance of 5 to 120 km downwind from instantaneous line-source releases. 

The maximum centerline concentration for the model in the Kth layer 
is given by the expression 

X CK { V y K = °’ Z K } = X K /{ LATERAL TERM) . (B-29) 


The average alongwind concentration is defined as 


*K °K /SjK 


(B-30) 


where the ground cloud passage time in seconds is 

Vk ~ 4 ‘ 3 (T xk/^K 


(B-31) 


The time mean alongwind concentration in the Kth layer is defined by the 
expression 


X K { V y K >Z K ' V 



(B-32) 
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where is the time in seconds over which concentration is to be averaged. 

The time mean alongwind concentration is equivalent to the average alongwind 
concentration when t equals T^. This complex set of relations, then, 

contains the computations performed in Model 3 to obtain the concentration- 
dosage mappings. 


48 



APPENDIX C 


INPUT PARAMETERS FOR THE NASA/MSFC 
MULTILAYER DIFFUSION MODEL 


There are two groups of input parameters for the model: the source 
input parameters which are vehicle and meteorologically dependent ( Table 
C-l) and the meteorological input parameters which are strictly dependent on 
meteorological conditions at launch time ( Table C-2) . These parameters 
include the special set employed in the layer breakdown model — Model 4. 

The source relationships given in Table C-l are determined in refer- 
ence to the stabilized ground cloud. The standard deviation of the crosswind 
source is 


cr 



Y(K) 

4.3 


(C-l) 


and the standard deviation of the alongwind source is 


cr 



MkI 

4.3 


(C-2) 


The source strength in the Kth layer is 


Q = 


f YlK) X{ K) \ 
v 2Y{k) X{K}/ 


Qn 


Z TK " Z BK 


(C-3) 


where Y{ K) and X{ K} describe the crosswind and alongwind dimensions of 
the cloud in the Kth layer, and is the total source strength in the ground 
cloud in units of mass. 

Equations (C-l) and (C-2) are based on the assumption that the along- 
wind and crosswind distribution of material in each layer is Gaussian and that 
the visible edge of the cloud represents the point at which the concentration is 
one-tenth the concentration at the cloud center in the Kth layer. Equation (C-3) 
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assumes that the cloud is spheroidal in the plane of the horizon and that the 
total source strength in the Kth layer is given by the relative cloud volume in 
the Kth layer. Because the models require the source strength per unit height, 
the total source strength in the Kth layer must be divided by the depth of the 
layer. 


The first nine meteorological parameters follow directly from the ther- 
modynamic and kinematic profiles of the atmosphere. The remaining two 
parameters (layer model) are empirical atmospheric constants. 


TABLE C-l. SOURCE INPUTS FOR THE MULTILAYER 
MODEL CALCULATIONS 


Parameter 


Layer 
Model: 1, 
2, 3 

Layer Break- 
down Model: 4 

1 

Definition 

Z R 

Z R 

Reference height in the surface layer 

Z BK 

Z BL 

Height of the layer base 

Z TK 


Height of the layer top 

t k 

t l 

Source (cloud) stabilization time 

X ryK 

X ryL 

Distance over which rectilinear lateral 
expansion occurs lownwind from an ideal 
point source 

o 

>> 

b 


Standard deviation of the crosswind source 
dimension in the Kth layer 

a {K> 
xo 


Standard deviation of the alongwind source 
dimension in'the Kth layer 


t* 

Time of layer breakdown 

Q 

^K 


Source strength in the layer 

J 


Scaling coefficient 
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TABLE C-2. LIST OF METEOROLOGICAL MODEL INPUTS 


Parameter 


Layer 
Model: 1, 

Layer Break- 


2, 3 

down Model: 4 

Definition 


U RL 

Mean wind speed at reference height 

*BK 

Z R 

Mean wind speed at the base of the 
layer 

^BL 

U TK 

U TL 

Mean wind speed at the top of the 
layer 

9 bk 

0 BL 

Mean wind direction at the base of 
the layer 

0 TK 

0 TL 

Mean wind direction at the top of the 
layer 

a AR^ T oK^ 

CT ARL^ T oL j 

Standard deviation of the wind azimuth 

angle at height z_, 

XV 

a ABK^ T oK^ 

a ABL^ T oL^ 

Standard deviation of the wind azimuth 
angle at the base of the layer 

^ATK^oK^ 

a ATL^ T oL^ 

Standard deviation of the wind azimuth 
angle at the top of the layer 

T oK 

T oL 

Reference time period 

°K 

“l 

Lateral diffusion coefficient 

TJ 

II 

p{L = l} 

Power-law exponent of the wind speed 
profile in the surface layer 


a ERL 

Standard deviation of the wind 
elevation angle at height z 

K 

Standard deviation of the wind 
elevation angle at the base of the Lth 
layer 

^EBL 



^ETL 

Standard deviation of the wind ele- 
vation angle at the top of the Lth layer 


^L 

Vertical diffusion coefficient 
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APPENDIX D 


TOXICITY CRITERIA 


A realistic evaluation of the potential hazards arising from high near- 
field concentrations of toxic effluents from solid rocket exhaust requires both 
a knowledge of the surface deposition of these effluents, which can be obtained 
with the NASA/MSFC Multilayer Diffusion Model (Appendix B), and toxicity 
criteria to evaluate the hazard from this surface deposition of effluent, which 
is the purpose of this discussion. The Federal Air Quality Criteria do not 
presently include any of the liquid or solid rocket exhaust effluents; however, 
the National Academy of Sciences does afford definite guidelines for the expo- 
sure to the toxic effluents associated with these exhausts. These guidelines 
are ecologically sound, based on the current limited knowledge of the effects 
of these effluents, and are the basis of the toxicity criteria that will be given 
[24, 25]. 

The primary effluents from any solid rocket exhaust are: aluminum 
oxide (A1 2 0 3 ), hydrogen chloride (HCl), carbon monoxide (CO), carbon 
dioxide (C0 2 ), hydrogen (H 2 ), nitrogen (N 2 ), and water vapor (Hp). While 
only the first four compounds are toxic in significant concentrations, there is 
always a potential hazard of suffocation from any gas which results in the 
reduction of the partial pressure of oxygen to a level below 135 mm Hg (18 
percent by volume at STP) . Oxygen level reduction does not appear to be a 
hazard from rocket exhaust because of the large volume of air that is entrained 
in these exhaust clouds; therefore, this potential hazard car be neglected in 
this discussion and attention can be directed to only the initial four toxic com- 
pounds. (A liquid rocket motor has only one toxic effluent — carbon monoxide.) 

The exposure level for toxic effluents is divided into three categories: 
public exposure level, emergency public exposure level, and occuptational 
exposure level. The public exposure levels are designed to prevent any detri- 
mental health effects both to all classes of human beings (children, men, 
women, the elderly, those of poor health, etc. ) and to all forms of biological 
life. The emergency level is designed as a limit in which some detrimental 
effects may occur, especially to biological life. The occupational level gives 
the maximum allowable concentration which a man in good health can tolerate — 
this level could be hazardous to various forms of biological life. 

The toxicity criteria for the toxic effluents in solid rocket exhausts are 
given in Table D-l. Public health levels for aluminum oxide are not g'ven 
because the experience with this particulate is so limited that, at best, the 
industrial limits are just good estimates. 


Hydrogen chloride is an irritant; therefore, the concentration criterion 
for an interval should not be exceeded [ 25] . Since hydrogen chloride is detri- 
mental to biological life, and in view of the fact that most launch sites are 
encompassed by wild life refuges, the emergency and industrial criteria for 
hydrogen chloride are not appropriate to the ecological constraints. Because 
of the large volume of air entrained in the exhaust cloud, the potential hazard 
from carbon monoxide and carbon dioxide can be, in general, neglected. 

Any detrimental health effects resulting from the combined toxicological 
action of these ingredients has been omitted because of a lack of knowledge in 
this area. However, investigations are currently underway to study this prob- 
lem and to learn more about the biological effects of hydrogen chloride. 
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TABLE D-l. AIR QUALITY TOXICITY STANDARDS' 
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